ANOMALY OF LINEARIZATION 
AND AUXILIARY INTEGRALS. 



BORIS KRUGLIKOV 

Abstract. In this note we discuss sonic formal properties of uni- 
versal linearization operator, relate this to brackets of non-linear 
differential operators and discuss application to the calculus of aux- 
iliary integrals, used in compatibility reductions of PDEs. 

Introduction 

Commutator [A, V] of linear differential operators A, V G Diff (vr, vr) 
in the context of non-linear operators F, G G diff (vr, vr) is up-graded to 
the higher Jacobi bracket {F,G}, which plays the same role in com- 
patibility investigations and symmetry calculusQ 

The linearization operator relates non-linear operators on a bundle 
vr with linear operators on the same bundle, whose coefficients should 
be however smooth functions on the space of infinite jets. The latter 
space is the algebra of ^-differential operators and we get the map 

£ : diff(vr,vr) ^ ^Diff(vr,vr) = C°°(J°°vr) ®c-(m) Diff(vr,vr), 

defined by the formula |KLV] 

ip(s)h = f^F{s + th)\t=o, F e diflF(vr,vr), s, /i G C°°(vr). 

However it does not respect the commutator: 

Example: Consider the scalar differential operators on M, so that 
vr = 1 and J°°(vr) = u,p = pi,p2, ■ ■ ■)■ Choose 

F = p^,G = p + c-x; {F,G} = 2cp =^ V,g} = 2cV,. 

If we commute ip = ^pV^ and ic = l^x, we get: [£^,^0] = —2p2T>x, 
so that we observe an anomaly. 

There are two reasons for this. The first is that the operator of 
linearization disregards non- homogeneous linear terms, which are im- 
portant for the Jacobi bracket. The second is the non-linearity itself. 
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The goal of this note is to discuss reasons and consequences of this 
anomaly (this also plays a significant role in investigation of coverings 
and non-local calculus |KKVj ). 

Acknowledgement. The results were obtained and systematized 
during the research stay in Max Planck Institute for Mathematics in 
the Sciences, Leipzig, in April-May 2007. 

1. Anomaly via Hessian 

The Jacobi bracket of non-linear operators F, G G diff(7r,7r) is ex- 
pressed via linearization as follows: 

{F, G} = ipG - icF. 

We also consider the evolutionary operators defined by duality: 

BpG = £gF. 

Since ic is a derivation in G, dp is a derivation (satisfies the Leibniz 
rule) and their union can be treated as the module of vector fields. 
These operators have no anomaly, i.e. the map 9 : C°°(J°^7r) 
Vect(J°°7r) is an anti-homomorphism: 

[9f, 9g] = -9{F,G}- 

This instantly implies Jacobi identity for the bracket {F,G}, so that 
(diff (vr, vr), {, }) is a Lie algebra [KLVj . 

The operators of universal linearization and evolutionary differenti- 
ation do not commute and this leads to the following 

Definition. The Hessian operator diff (vr, vr) x diff (vr, vr) ^ Diff (vr, it) 
is defined by the formula 

ResspG=[9G,iF]. 

We will also write ResspiG, H) = Ressp G{H) for F,G,H e diS{n, n) 
and note that Hessi? = for linear operators F, because in this case 
ip = F, which reduces the claim to the commutation of left and right 
multiplications. 

Next we note that the Hessian HesSi? is symmetric: 

Lemma 1. HesspiG, H) = Hessi7(if, G). 
Indeed: 

HessF(G, H) = dcipH - £f9gH = 9g9hF - If^hG, 

so that 

Hessir(G', H) - ResspiH, G) = [9g, 9h]F - If{H, G} 

= -d{G,H}F-£F{H,G} = Q. 
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Now we can express the anomaly of linearization via the Hessian: 
Proposition 2. [£p, £g] — i{F,G} = Hessc F — Hess^ G. 
Indeed we have: 

[epjaW^eFdHG-eGdHF 

= Bh('^fG - ioF) - Ress f(H, G) + HessG(i/, F) 
= dH{F, G} - HessF(G', H) + HessG(F, H) 

= lt{F,G}H + (HessG F - Hess J. G)H. 

Finally let us express the Leibniz identity for non-linear operators 
and the Jacobi bracket. For linear operators it is well-known, but for 
non-linear ones there's an anomaly: 

Proposition 3. {F, IqH} = 1{f,g}H + Ig{F, H} - HessF(G', H). 
This is obtained as follows: 

{F,eGH} = eFiGH-9FiGH 

= [iF, iG]H + iGi^F - 9f)H - Hessc(F, H) 

= hF,G}H + tG{F, H] - Hess^(G, H). 

2. Coordinate expressions 

A local coordinate system (a;\ u^) on tt induces the canonical coordi- 
nates {x\pi^) on the space J°^7r, where a = {ii, .... i„) is a multi-index 
of length \cr\ = ii + ■ ■ ■ + in- The operator of total derivative of multi- 
order a (and order \a\) is = T>^^ ■ ■ ■ T>1^, where "Dj = d^-i+Y^ p-^^^^d^j . 

The linearization of F = (Fi, . . . , Fr) is If = (K^i), ^{Fr)) with 

where V^^ denotes the operator V„ apphed to the j-th component of 
the section from C°°(7r). 

The i-th component of the evolutionary differentiation do corre- 
sponding to G = (Gi, . . . , Gn) equals 

where dj denotes the operator applied to the i-th component of 

Pa Pa 

the section from C°°(7r). 
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Then i-th components of the Jacobi bracket is given by 

{F,Gh = ■ - ^'^(^^■) ■ 

These formulas are known |KLV] . It is instructive to demonstrate the 
Jacobi identity in coordinates. For this we need the following assertion. 

Lemma 4. In canonical coordinates on J'^tt: 

(the difference of multi-indices a — x is defined whenever x d a), 
the summation is by x counted with multiplicity. More generally for 
vector differential operators ifV^J^ is the operator!)^ acting on thej-th 
component, then the above formula holds true for such specification. 

This follows from iteration of the formula \d 3 , VA =83 . Thus 

{F,{G,H}} = Y,F,V,-.{G,^)V,+.{H)-F,V^_^{H,^)V^+^{G) 
- G,^,MH)V^iF) + H,^,MG)V^{F) 

- {Gp„V^^^{Hp^_J - Hp^V^^^{Gp^_J)Vr{F), 

which yields Ecyciic{^, {G, H}} = 0. 
Now we write the Hessian: 

Hess^(G', H) = J2 Fp^PrT^aG ■ V^H, 

and its symmetry in G, H and vanishing for linear F is obvious. 
The compensated Leibniz formula can be written as follows: 

{F,IgH}-1{f,g}H~Ig{F.H} = 
Y,Fp.V^-.{Gp^)Vr+^{H)-{Gp^p^V,{H)V,{F)+Gp^dp^V,{H))V„{F) 
-{F,^,MG)+Fp^dpMG))VAH)HGp.pMF)+Gp^dpMF))VAH) 
-Gp^ {V,^^{Fp^)Vr+^{H)-V,,^{Hp^)Vr+^{F)) = - Hessi.(G, if) 
and the anomaly in commuting linearizations is: 

[^f,^g]- W} = 

Fp^V^^^iGp^) V^UH) - Gp^V^^^iFp^) V^^H) 
-{Fp^pMG)+Fp^dpMG))Vr{H)+iGp^pMF)+Gp^dpMF))Vr{H) 

= HessG(F, H) - HessF(G, H). 
This gives an alternative proof of Propositions [3] and [21 



3. Auxiliary integrals 



Definition. An operator G G diff(7r,7r) is called an auxiliary integral 
for F e diff(7r,7r) if 

{F,G} = hF + i^G 

for some operators X G diff(7r,7r) and fi ^ '^Diff (tt, tt) ■ F \ {0}. The 
set of such G is denoted by Aux(F). 

It is better to denote Aux^(F) the space of G satisfying the above 
formula with some fixed /i G diff(7r, vr), because it is a vector space. 
Then Aux(F) = U^Aux^(F). We can assume ord(/i) < ord(F) for 
scalar operators, i.e. rankvr = 1. 

With certain non-degeneracy condition for the symbols of F, G the 
following statement holds: 

Theorem 5. A non-linear differential operator G is an auxiliary inte- 
gral for another operator F iff' the system F = 0,G = is compatible 
(formally integrable). 

The generic position condition for the symbols of F, G is essential. If 
vr = 1 is the trivial one-dimensional bundle, this condition is just the 
transversality of the characteristic varieties Char'^(F) and Char'''(G') 
in the bundle P'^T*M (after pull-back to the joint system F = G = 
in jets); in this form it is a particular form of the statement proved 
in IKL2I . For rankvr > 1 the condition is more delicate and will be 
presented elsewhere. 

Notice that Auxo(-F) = Sym(F) is the space of symmetries of F. 
This is a Lie algebra with respect to the Jacobi bracket. It can be 
represented as a union of spaces 

Sym,(F) = {H : ipH = ie+nF}, 9 G difT(7r, tt), 

which are modules over SymQ(F). More generally we have the graded 
group: Symg,(F) + Symg„(F) C Symg,+g„(F) 
Let us assume G G Aux^(F), H G Symg(F), i.e. 

{F, G] = £xF + i^G, {F, H} = (eF 

Then denoting adj^ = {H, ■} = Eh — 9h we get: 
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adpiG, H] = {adp G, H} + {G, adj. 

= -{H4xF + ^^G] + {G,^eF] 

= W}^+^A{i^, i/}+HessH(A, F)+^{^,;,}G+£^{G, i/}+Hess^(/., G) 

- ^{^^,G}i^ - ^e{F, G] - HessG(^, F) 

= {k\H} + [^A, ^e] - k6,G} + Hessj^ A - HessG e)F + f^G, H] 

+ - ^e^M + Hessi/ /x)G. 

Thus {G, H} is an auxiliary integral for F if = i{fj,,H} + Hessj^ /i 
(the "iff" condition means the difference annihilates G), which can be 
written as 

II e Ker[(£e + - Hessj^) o £]. 

Such a pair 9 e sym*(F) = Sym(F)/ Symo(F), H e Symg(F) deter- 
mines the action of the second component 

adH : Aux^(F) Aux^(F). 

Also since 

- 9{G,H}(J' - Hess//(G, /x) = - adn dailA - fiess^lA*, G) - £^{G, H}, 
we have: 

£^{G, H] + (£|^„H} - + Hess^ /i)G = -(ad^ +4)^^G. 

Thus ilH e Symg(F), i.e. (ad// +£e)F = 0, and e Ker[(ad// +£e)o-^], 
i.e. (ad/j- -l-^e)^// = 0, then 

ad// : Aux^(F) ^ Sym(F). 

4. Symmetries and compatibility 

It has been a common behef that if G e Sym(F), then the system 
F = 0, G = is compatible, which forms the base of investigation for 
auto-model solutions. This is however not always true. 

Example: Let F, G be two linear diagonal operators with constant 
coefficients. Then {F, G} = (in this case the Jacobi bracket is the 
standard commutator), so that G is a symmetry of F. However the 
system F = 0, G = is usually incompatible: for generic F, G of 
the considered type the only solution will be the trivial zero vector- 
function. 

More complicated non-diagonal operators are possible, but it would 
be better to consider non-homogeneous linear operators. Then if the 
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coefficients are constant and generic, the linear matrix part commute, 
but the system F = 0,G = may have no solutions at all. 
For instance if we take 
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then {F, G} = 0, so that G G Sym(F), while the system F = 0,G 
is not compatible, and moreover its solutions space is empty. 

Thus the flow Ut = G{u) on the equation F = has no flxed points 
(no auto- model solutions). Here t is an additional variable (x is the 
base multi- variable for PDEs F = and G = 0), so that G G Sym(F) 
can be expressed as compatibility of the system 



F{u) = 0, ut = G{u) 



0, 



I.e. 



while symmetric solutions correspond to the stationary case Ut 
compatibility of the system F{u) = 0, G{u) = [0. 

However if the non- degeneracy condition assumed in Theorem [5] is 
satisfled, then auto-model (or invariant) solutions exist in abundance, 
namely they have the required functional dimension and rank as Hilbert 
polynomial (or Cartan test [C]) predicts, see |KL. 



J4 . 



Remark. Symmetric solutions are the stationary points of the evolu- 
tionary fields and they are similar to the fixed points of smooth vector 
fields on M", which must exist provided the vector field is Morse at 
infinity. The non- degeneracy condition plays a similar role. 

Many examples of auto-model solutions and their generalizations can 
be found in |BK[[OlllOvj . non-local analogs use the same technique and 
similar theory pOM [KKllKKVj . 

Compatible systems correspond to reductions of PDEs and are some- 
times called conditional symmetries by analogy with finite-dimensional 
integrable systems on one isoenergetic surface [FZj. But the rigorous 
result must rely on certain general position property for the symbol 
of differential operators, otherwise it can turn wrong IKL2I [KL3] . The 
method based on this approach makes specification of the general idea 
of differential constraint and is described in |KLi| . 



I am grateful to S.Igonin and A.Verbovetsky for an enlightening discussion 
about the results of |KL2[ IKL3I and the symmetry condition. 
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5. Conclusion 



In this note we described the higher-jets calculus corresponding to 
symmetries and compatible constraints, basing on the Jacobi brackets. 
Another approach to integrability of vector systems is given by minimal 
overdetermination and it uses multi-brackets of differential operators 

{■■■}: A™+i difr(m ■ 1, 1) ^ difr(m ■ 1, 1) 

introduced in IKL3I , which are governed by the non-commutative Pliicker 
identity. 

Following this approach a minimal generalization of symmetry for 
F = (Fi, . . .,Fm) e difr(7r,7r) with vr = m - 1 is such G G diff(7r, 1) that 

{Fi, . . . , G} = ie,Fi + ■■■ + ie^F^. 

With certain non-degeneracy assumption IKL3I this implies that the 
overdetermined system F = 0, G = is compatible (formally inte- 
grable). 

A more advanced algebraic technique would yield another higher- 
jets calculus producing anomaly that manifests in non- vanishing of the 
expression 

Implications for vector auxiliary integrals and generalized Lagrange- 
Charpit method follow the same scheme. 

References 

[BK] G. W. Bluman, S. Kumei, Symmetries and differential equations, Appl. 
Math. Sci. 81, Springer, 1989. 

[C] E. Cartan, Les systemes differentiels exterieurs et leurs applications 
geometriques (French), Actuahtes Sci. Ind. 994, Hermann, Paris (1945). 

[FZ] W.I. Fushchych, R.Z. Zhdanov, Conditional symmetry and reduction of par- 
tial differential equations, Ukrain. Math. J. 44 (1992), 970-982. 

[KKV] P. Kersten, I. S. Krasilschik, A. Verbovetsky, Hamiltonian operators and 
£* -coverings, J. Geom. and Phys., 50 (2004) 273-302. 

[KLV] I. S. Krasilschik, V. V. Lychagin, A.M. Vinogradov, Geometry of jet spaces 
and differential equations, Gordon and Breach (1986). 

[KLi] B. S. KrugHkov, V. V. Lychagin, A compatibility criterion for systems of 
PDEs and generalized Lagrange- Charpit method, A. LP. Conference Proceed- 
ings, Global Analysis and Applied Mathematics: International Workshop on 
Global Analysis, 729, no. 1 (2004), 39-53. 

[KL2] B. S. Krughkov, V. V. Lychagin, Mayer brackets and solvability of PDEs - 
II, Trans. Amer. Math. Soc. 358, no.3 (2005), 1077-1103. 

[KL3] B. S. KrugHkov, V. V. Lychagin, Compatibility, multi-brackets and inte- 
grability of systems of PDEs, prepr. Univ. Troms0 2006-49; ArXive: 
[math.DG/0610930 

8 



[KL4] B. S. Kruglikov, V. V. Lychagin, Geometry of Differential equations, in: 
D. Krupka, D. Saunders, Handbook of Global Analysis (2007); prepr. 
IHES/M/07/04. 

[KK] I.S. Krasilshchik, P.H.M. Kersten, Symmetries and recursion operators for 
classical and super symmetric differential equations, Kluwer (2000). 

[01] P. Olver, Applications of Lie groups to differential equations, Graduate Texts 
in Mathematics, 107, Springer- Verlag, New York (1986). 

[Ov] L. V. Ovsiannikov, Group analysis of differential equations, Russian: Nauka, 
Moscow (1978); Engl, transl.: Academic Press, New York (1982). 

Institute of Mathematics and Statistics, University of Troms0, Tromso 
90-37, Norway. 

E-mail address: kruglikov@math.uit.no 



9 



